Nonlinear subdiffusive fractional equations and aggregation phenomenon 
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In this article we address the problem of the nonlinear interaction of subdiffusive particles. We 
introduce the random walk model in which statistical characteristics of a random walker such as 
escape rate and jump distribution depend on the mean field density of particles. We derive a set 
of nonlinear subdiffusive fractional master equations and consider their diffusion approximations. 
We show that these equations describe the transition from an intermediate subdiffusive regime to 
asymptotically normal advection-diffusion transport regime. This transition is governed by nonlinear 
tempering parameter that generalizes the standard linear tempering. We illustrate the general results 
through the use of the examples from cell and population biology. We find that a nonuniform 
anomalous exponent has a strong influence on the aggregation phenomenon. 

PACS numbers: 



I. INTRODUCTION 

Anomalous subdiffusion is a widespread phenomenon 
in physics and biology [1H3J. It is observed in the trans- 
port of proteins and lipids on cell membranes [4] , RNA 
molecules in the cells Q, signaling molecules in spiny 
dendrites |6[, and elsewhere. Apart from fractional Brow- 
nian motion, the linear fractional equations are the stan- 
dard models for the description of anomalous subdiffusive 
transport Q- In these models the diffusing particles do 
not interact. The question then arises as to how to extend 
these equations for the nonlinear case, involving particles 
interactions. These nonlinear effects are typical and very 
important in cellular and population biology. It is well 
known that many biophysical processes in microorgan- 
isms depend on the their population density. A typical 
example is the quorum sensing phenomenon, in which 
microorganisms coordinate their behaviour according to 
their local population density [7]. Other examples are 
a cellular adhesion, which involves the interaction be- 
tween neighbouring cells [8Hll|: and the volume-filling 
effect, which describes the dependence of cell motility on 
the availability of space in a crowded environment [12j . 
The understanding of macroscopic phenomena like cell 
and microorganizm aggregation requires an understand- 
ing of how individual species interact through attractive 
or repulsive forces (see, for example, [13| and references 
therein). The attraction between individuals may result 
from various social interactions such as mating, settle- 
ment, defense against predators, etc. While the repul- 
sion ma y o ccur due to low resources in highly populated 
regions [13, [lj] . Note that microorganisms interact both 
directly and indirectly via signaling molecules. 

The main purpose of this paper is to incorporate these 
nonlinear effects into subdiffusive equations. Our aim 
is to take into account the interaction between particles 
on the mesoscopic level, at which the random walker's 
characteristics depend on the mean field density of par- 
ticles. Our intention is to derive the nonlinear subd- 
iffusive fractional equations for the density of particles 
and apply these equations to the problem of aggrega- 
tion. In this paper we use two different approaches that 



are based on the density-dependent dispersal kernels and 
density-dependent jump rate. Note that several theoreti- 
cal studies have been devoted to nonlinear generalizations 
of linear subdiffusive fractional equations. However, most 
research has been focused on the problem how to incor- 
porate the nonlinear reactions into subdiffusive equations 
|L5l - |24] . The aim of this paper is to study the nonlinear 
subdiffusive transport processes involving the anomalous 
trapping of particles and their interactions. 

In this paper we deal with the random walk model in- 
volving a residence time dependent escape rate and the 
structural density of particles. This has been used by 
many authors for the analysis of the non-Markovian ran- 
dom walks @, Qji [la, UM, [23 ■ It turns out that this lin- 
ear model is the most suitable for further nonlinear gen- 
eralizations. We consider a 'space-jump' random walk 
in one space dimension. The particle waits for a ran- 
dom time (residence time) T x at point x in space before 
making a jump to another point. The random residence 
time T x is determined by the probability density function 
ip(x, t) — Pr {r < T x < t + dr} . The key characteristic 
of this random walk is the escape rate 7 from the point 
x. It depends on the residence time r and the position 
x : 7 = 7(2;, r). This rate can be rewritten in terms of 
the probability density function i/j(x,r) and the survival 
probability ^/(x, r) = J t °° ip(x, u)du as follows |27j 



l(x,r) 



^(x,t) 
*(x,r)' 



(1) 



Let £(x, t, t) be the structural density of particles at 
point x at time t whose residence time lies in the interval 
(t, t + dr) . We assume that at t = all particles have 
zero residence time: £(x,T,0) = po{x)S(r). The mean 
escape rate, i (x, t) , of particles from the point x can be 
defined as 



i(x,t) 



•y(x,T)£(x,T,t)d,T 



and the density of particles p {x, t) 

p(x,t) = / £(x,T,t)dT. 

Jo 



(2) 



(3) 



One of the main results in the non-Markovian random 
walk theory is that the mean escape rate i (x, t) can be 
written as a convolution [3j 



i (x,t) 



K (x,t — t) p (x, t) dr. 



(4) 



Here K[x, t) is the memory kernel defined by its Laplace 
transform 



K(x,s) 



4> (x,s)" 



(5) 



where tp (x, s) 



r- 



e st tp(x, t)cIt and * (a;, s) 



L e~ st fy(x,T)dT [3|. The master equation for the den- 
sity p (x, t) can be written as the balance of particles at 
the point x 

— = / i (x — z, i)w (z\x — z) dz — i (x, t) , (6) 

vt Jr 

where w (z\x) is the dispersal kernel. 

In the anomalous subdiffusive case, the survival prob- 
ability Vfr.j, (x, t) can be modelled by the Mittag-Lcfficr 
function 128] 



V, i (x,t) = E ll 



tq(x) 



fj,(x) 



, < u(x) < 1, (7) 



spatial variations of fractal exponent u {x). This leads to 
the anomalous aggregation of particles [26| . In heteroge- 
neous biological systems, in which the exponent p(x) is 
space-dependent, the question arises as to whether this 
anomalous aggregation of a population can be prevented. 
Therefore it is an important problem to find the way how 
to regularize subdiffusive fractional equations. One way 
is to incorporate random killing, which ensures regular 
behaviour in the long-time limit [30J|. The aim of this 
paper is to address this problem through a nonlinear es- 
cape rate that takes into account repulsive forces between 
particles. In the next section we derive the nonlinear gen- 
eralization of the fractional equations like ([9]) and (TTOl) . 



II. NONLINEAR SUBDIFFUSIVE 
FRACTIONAL EQUATIONS 

There are two major ways in which nonlinear density- 
dependence effects can be implemented into non- 
Markovian and subdiffusive transport equations. The 
simplest way is to take into account the dependence of 
jump density w on p. This dependence can take into ac- 
count various non-linear effects such as adhesion, quorum 
sensing, volume filling, etc. However, we begin with more 
complicated case of the random walk model for which the 
escape rate is a function of the residence time and the lo- 
cal density of particles. 



where p(x) is the space-dependent anomalous exponent. 
In what follows we assume for simplicity that To is con- 
stant. The Laplace transform of the memory kernel 
Kf,(x,t) is K^{x,s) = s i-M*) To -M*) an d the integral 
anomalous escape rate i (x,t) can be written as 



i {x,i) = 



1 



T() 



fj.(x) 



■Dl~ Kx) p(x,t) 



(8) 



l — n(x) 

where T> t is the The Riemann-Liouville derivative 

with varying anomalous exponent [29|. Substitution of 
([8]) into ((6]) gives the integral fractional equation with 
space-dependent anomalous exponent 



dp 

at 



i 



-V 



l — fj,(x — z) 



to 



TQt j.(x-z) 
p,(x) t 



p{x — z,t)w (z\x — z) dz 
p{x,t). (9) 



Using the Taylor series expansion in terms of z, and a 
symmetric dispersal kernel w, we obtain the subdiffusive 
fractional equation 



dp 
dt 



dx 2 



D^(x)V 



l-H(x) 



p(x,t) 



(10) 



with the fractional diffusion coefficient D^x) — 
a 2 (x)/T ^ x \ where a 2 (x) — J z 2 w{z\x)dz. 

It has been found recently that subdiffusive fractional 
equations are not structurally stable with respect to the 



A. Nonlinear escape rate with memory and 
nonlinear tempering 

The main problem with this anomalous escape rate 
([8]) is the phenomenon of anomalous aggregation [25|, [26j . 
Nonuniform distribution of the anomalous exponent p(x) 
over the finite domain [0,L] leads to 



p(x,t) — > 5{x — Xm) as 



t — > 00. 



(11) 



Here xm is the point in space where the anomalous ex- 
ponent p (x) has a minimum. The problem is that that 
the escape rate i (x,t) is a linear functional of the den- 
sity of particles p(x,t) and does not take into account 
nonlinear effects of repulsive forces which, in many sit- 
uations, can prevent anomalous aggregation. According 
to (jTT|) all particles aggregate into a small region around 
the point x — xm forming a high density system. To pre- 
vent such anomalous aggregation, one can assume that 
the overcrowding leads to an increase of repulsive forces 
and a corresponding correction of the anomalous escape 
rate 7(3;, r). 

We assume that the probability of escape due to the 
repulsive forces is independent from anomalous trapping. 
We define this probability for a small time interval At as 



a(p(x,t))At + o(At). 



(121 



Here a(p) is the transition rate which is an increasing 
function of the particles density p. Another interpreta- 



tion can be given in terms of the quorum sensing phe- 
nomenon [3] . The large cell density can lead to the local 
over-depletion of nutrients and oxygen and as a result 
cells can change phenotype from a proliferating state to 
a migrating one 21]. In another words, when the concen- 
tration of cells is low, a(p) = 0, but if the concentration 
of cells p(x,t) reaches a certain level p cr : a(p) ^ for 
P > Per- Of course, one can assume a non-monotonic de- 
pendence of the transition rate a(p) on p. For example, 
at low cell densities a(p) could decrease with p, while 
at high densities a{p) could be an increasing function: 
a(p) = a (l - aip + a 2 p 2 ) [H, [l| . 

Firstly we formulate the Markovian model for the ran- 
dom walk of particles with the density-dependent escape 
rate. The balance equation for the structural density 
£(x, r, t) takes the form 

£(x,r + At,t + At,) = £(x,T,t){l-j(x,T)At) x 

(l-a(p(x,t))At)+o(At). 

where 1 — "f(x, r)At is the survival probability during 
time At due to the trapping, and 1 — a(p(x,t))At is the 
survival probability during time At corresponding to the 
repulsion forces between particles. In the limit At —> 0, 
we obtain 



dt dr 



[l(x,r) + a{p{x,t))]£, 



(13) 



where the effective transition rate is the sum of two es- 
cape rates: 



7(z,t) + a(p(x,t)). 



(14) 



The second term a(p(x, t)) can be treated as the cor- 
rection of the basic escape rate 7(2;, r) due to repulsive 
forces. The boundary condition for £(x,t, t) at zero res- 
idence time t = is 



£(x,0,t) = / / [j(x-z,T) + a(p(x-z,t))] x 

£,(x — z, t, t)w (z\p(x — z, t)) drdz. (15) 
To obtain the effective escape rate 

i(x,t)= / j(x,T)£(x,T,t)d,T + a(p)p(x,t) (16) 



we solve lfT3")) by the methods of characteristics and sub- 
stitute £,(x,t, t) into (fl~5|) and p^|) . Finally we obtain 



i(x,t) = / K(x,t-T)e-f> Mx ' s))ds p{x,T)dT 



+a(p)p(x,t). 



(17) 



One can see that the first term in (fT7f involves the expo- 
nential factor, e~ ^ <x(p(x,s))ds ^ wn j cn can }-,£ interpreted 
as a nonlinear tempering. The main feature of this mod- 
ified escape rates is that although the local rates 7(0;, r) 
and a(p(x,t)) are additive (see (fT4l ). the corresponding 



terms in the integral escape rate (jTTJ) are not additive. 
This is clearly non-Markovian effect. One of the main 
aims of this paper is to find out what the implications 
of this effect on the long-time behaviour of the density 
p(x,t) are (see section III). 

Note that when 7 does not depend on the residence 
time variable r, the survival function has an exponential 
form ^(x, r) = e -7 ^ 1 ", the Laplace transform K (x. s) = 
j(x) and the Markovian escape rate 



i(x,t) = [7(2:) +a{p)]p(x,t). 



(18) 



Substitution of this into © gives the closed nonlin- 
ear Kolmogorov- Feller equation for p(x,t) [3>]. One can 
find various nonlinear diffusion approximations of this 
equation assuming density-dependent dispersal kernels 
w (z\p (x, t)) and density-dependent jump rate a(p (x, t)) 



B. Anomalous subdiffusive case 

In the subdiffusive case, it follows from (IT71) that in- 
stead of (|8j) we have 



i {x,t) 



+a(p)p(x,t), 



M*,t) 



p(x,t) 



where 



$(x,t) 



a (p(x, s)) ds. 



(19) 



(20) 



This function plays a very important role in what follows. 
For the jump density involving only two outcomes 

w (z\x. t) = r(x, t)8(z - a) + l(x, t)8(z + a), (21) 

where a is the jump size, we obtain the master equation 
for p (x, t) in terms of the total flux i (2, t) 



dp , ... . 

— = r{x — a,t)i(x ~ a 7 t) 

+l(x + a, t)i{x + a,t) — i(x, t). 



(22) 



Here r(x,t) is the jump probability from the point x to 
x + a, l{x, t) is the jump probability from the point x to 
x — a. One can introduce the potential U(x, t) that induce 
the bias of the random walk in one direction [3l|, |32[ . In 
the case of chemotaxis, U{x, t) can be interpreted as the 
density of chemotactic substance [33, [34J. In the limit 
a — > 0, by using 



dU 

r(x — a,t) — l(x + a,t) = —(3a— h o(a) 7 

ox 



we obtain 



dp = _, 3a 2JL 
dt dx 



dU . ' 
—i{x,t) 
ox 



^ 2 d 2 i(x,t) 

dx 2 



(23) 



(24) 



This equation can be rewritten in the form of a nonlinear 
fractional Fokker-Planck equation 



dp = o a 2_^_ 
dt dx 



d 2 



dU ( e~* 



^( i ^MA'-*»[eM + «(* 



9a; 2 



270^) 



2^-/*(-)[e*p] + (p)p 



(25) 



This equation is one of the main results of this paper. 
It describes the transition from subdiffusive transport 
to asymptotic normal advection-diffusion transport. We 
discuss this asymptotic transition in section III. Note 
that the above equation can be rewritten in terms of the 
fractional diffusion coefficient 



DM 



which is finite in the limit a — > and To 



(26) 



0. 



C. Nonlinear jump distributions 

In this subsection we discuss various nonlinear jump 
distributions leading to nonlinear fractional equations. 
As we mentioned previously, this is the simplest way 
to incorporate nonlinearity into subdiffusive fractional 
equations [26]. As long as the escape rate i (x,t) is de- 
termined, we can define "where to jump" through the 
dispersal kernel. Let us assume that the dispersal kernel 
w depends on the density of particles: w — w(z\p(x,t)). 
In what follows we restrict ourselves to the subdiffusive 
case for which the rate i(x,t) is determined by (jSJ. 



One can also introduce the nonlinear drift term generated 
by Bernoulli jump distribution [131 ] 



w(z\p) = -(l + au(p))8(z-a) + 
~(l-au(p))5(z + a), 

where the positive function u (p) takes into account the 
fact that the jump on the right is more likely than a jump 
to the left. It is assumed that 1 — au (p) is not negative; 
that is, u (p) < a~ 1 . The fractional master equation takes 
the form 

9p_ = l + au{p(x-a,t)) vl (x _ a) 

dt 2T v{x-a) * PV ; 

l-Ptt(p(s + a,t)) i_ M(a+o) , 

2r ^+°) * P{X + a,t) 



TQ 



l —V\-^p{x,t) 



fi(x) 



(28) 



Using the Taylor series expansion in (|28p. we obtain the 
fractional Fokker-Planck equation 



dp d 

dt dx 

dx 2 



2D^x)u{p{x,t))V] Kx) p 
D,{X)V]-» {X) P 



(29) 



with the nonlinear advection term and fractional diffu- 
sion coefficient D^{x). 



2. Nonlocal nonlinearity 



1. Local nonlinearity 

First we consider the case for which the jump depen- 
dence on p is local. For example, one can use Gaussian 
dispersal kernel with rapidly decaying tails as 



w{z\p) 



1 



V^Hp) 



cxp 



2a 2 (p) 



An increasing dispersion a 2 (p) describes the effect of 
the local repulsive forces, while a decaying c 2 (p) cor- 
responds to attractive forces. Using the Taylor series 
expansion in terms of z in the master equation (|S|) with 
w (z\p(x — z, t)), we obtain the subdiffusive equation 



dp 

dt 



dx 2 



,i-^0) 



K,( P )vr^>p(x,t) 

with the nonlinear fractional diffusion coefficient 

°\P) 



(27) 



d„(p) 



TQ 



fi(x) ' 



Here we consider the volume filling effect and cell-to- 
cell adhesion [12[ . One can use the jump density 



w (z\p) = r(p)5 (z - a) + l(p)5 (z + a) 



(30) 



where r(p) is the probability of jumping from x to x + a 
and l(p) is the probability of jumping from x to x — a 
defined as follows 



r(p) 



HP) 



X(p) + p{p) 



, l{x\p) 



m(p) 



Kp) + Kp) ' 



(31) 



Here 



p(p) = ^^-IPi^+a.t) e - m p{x-a,t) ^ ( 32 ) 



where the factor e - m p( x + a t) [ s U ged to model volume fill- 
ing: the jump probability to the right r(p) is reduced by 
the the presence of cells at the point x + a; m is the vol- 
ume filling parameter. The factor e~ qptyX ~ a ^ describes 
the adhesion effect, which says that the jump probability 
to the right r(p) is reduced by the presence of cells at 



the point x — a; q is the adhesion coefficient. The master 
equation for the mean density takes the form 



.©l-foo^f) 



(33) 



III. TRANSITION TO ASYMPTOTIC 

ADVECTION-DIFFUSION TRANSPORT 

REGIME AND AGGREGATION PHENOMENON 

In this section we discuss a transition from a subd- 
iffusivc transport regime to an asymptotically normal 
advection-diffusion transport regime. This transition is 
governed by the nonlinear fractional equation (125[) . It 
involves the exponential factor e~*' :c, *' ) , with $ (#, t) = 
J a(p(x,s)) ds, that can be considered as a nonlinear 
tempering. It generalizes the linear tempering, where 
the power law waiting time distribution is truncated by 
an exponential factor exp(— at) [351 ] . We expect that in 
the limit "J? (x, t) —> oo we recover the classical Fokker- 
Planck equation, while in the intermediate asymptotic 
regime, when $ (x, t) << 1 we have a subdiffusive trans- 
port. If we take the limit $ (x,t) — > in (|25p for small 
a, we obtain the standard fractional Fokker-Planck equa- 
tion ED, El 



dp 


d 


dt 


dx 


+ 


d 2 
dx 2 



2D,(x)f3^-Vl-^p 
D^X)V\-» {X) 



(34) 



where D^(x) is defined by (|26[) . Nonuniform distribu- 
tion of the anomalous exponent p(x) in (|3"4"]l leads to 
the aggregation of particles at the point of minimum of 
//(.t) [19]. This phenomenon has been observed in exper- 
iments on phagotrophic protists when "cells become im- 
mobile in attractive patches, which will then eventually 
trap all cells" [33]. In this case the potential field U(x) is 
irrelevant to the distribution of p, and anomalous expo- 
nent p{x) dominates. However, this aggregation around 
one point is just a transient phenomenon. If we consider 
the long time limit t — ¥ oo then we obtain the classical 
Fokker-Planck equation (see next two subsections). The 
stationary solution of this equation represents ultimate 
cells aggregation determined by the potential field U(x) 
and the spatial distribution of anomalous exponent p{x). 



A. Aggregation in the linear case 

Firstly, we assume that the local escape rate a is inde- 
pendent of time, such that <J> (a;, t) = a(x)t. The Laplace 



transform of i (x, t) , defined by (fT9]l . gives 

(s + a(x)) l ~ fl( - x ^ 

i(x,s) = ^-j p(x,s)+a(x)p(x,s). (35) 



In the limit s — > (t — > oo) the dependence of the first 
term in (|35[) on s disappears, so the total escape rate 
i (x, t) can be written in the Markovian form 

i(x,t) = X fi (x)p(x,t) , 

where the effective rate of escape A M (a;) is 

a(x) 



<M X ) 



(T a(x)) 



li(x) 



a(x). 



The essential feature of this rate parameter A M (x) is that 
it depends on the fractal exponent p{x). This escape rate, 
together with (f24|) . leads to the standard Fokker-Planck 
equation 



dp 

lit 



dU 
2f3—D(x)p(x,t) 



+l&[D(x) P (x,t)} 



d_ 

dx 

d 2 

dx- 



(36) 



where the diffusion coefficient D (x) depends on p(x) and 

a(x) : 



D(x) = 



a 2 A p (x) 



(T a(x)) 



l-fi(x) 



T a(x) 



2 2r 

When the product r^a is small, the term (tqOc) 
dominant in the anomalous case p[x) < 1, so Z? M (x) can 
be approximated as 



is 



D(x) 



a 2 (roa(x)) 1 -^ 



2r 



Let us find the stationary solution p st (x) to ([36]) in the 
interval [0,L]. We use the reflective boundary conditions 
at x = and x = L which guarantees the conservation 
of the total population. We introduce the new function 

p(x) = D(x)p st (x). 

It follows from (l36t that this function obeys the equation 

d 2 p(x) 



d_ 
dx 



2p g P{x) 



dx 2 







(37) 



with the solution in the form of the Boltzmann distribu- 
tion 

p(x) = N- 1 cxp [-2fiU(x)\ . 

Thus, the steady profile is 

p st {x) = N^D- 1 (x) exp [-20U(x)] , (38) 




FIG. 1: Stationary profiles p s t(x) for the linear distribution 
of the potential U(x) — 5x (solid line) and U(x) — (dashed 
line), = 1, no = 0.9, a = f0" 4 , k = 2.f9. 



where N is determined by the normalization condition 
N = J L D- 1 (x) exp [-20U(x)] dx. 

To illustrate how the nonhomogeneous anomalous ex- 
ponent p(x) affects the aggregation pattern we consider 
the steady profile (|38|) in the interval [0, 1] for two cases: 
(a) U(x) = and (b) U{x) = mx for which the exponent 
p(x) has the form 

p{x) — hq exp (—kx) , 

where /iq < 1 . Fig. I shows that for the uniform potential 
U(x) — 0, the particles have the tendency to aggregate in 
the region of the small values of p(x) (dashed line). When 
the potential field forces the particles to move from the 
right to the left (U(x) = 5a;), the steady profile is not 
monotonic (solid line). 



B. Aggregation in the nonlinear case 

Now we consider the nonlinear case when the escape 
rate a depends on the density p. For living systems, this 
nonlinear dependence results from a coupling between 
the local density of cells and the intensity of the response 
of individual cell to external signals. We assume that in 
the limit t —y oo the stationary distribution p s t(x) exists. 
Then as t — » oo, the tempering factor e - *^'*) can be 
approximated by 

e -oc{p,t(x))t _ 

The stationary escape rate i s t(x) corresponding to (jTTJ) 
can be written in terms of the Laplace transform K (x, s) 
as follows 



i st (x) = K(x,a(p st )) + a(p at ) Pst(x) 



(39) 



This steady rate i s t (x) has the Markovian form in which 
the rate parameter consists of two terms K (x,a(p s t)) 
and a(p s t). The dependence of the first term on a(p s t) is 
due to the non-Markovian character of transport process. 
This effect does not exist in the Markovian case for which 
if is a function of x only. 



For the subdiffusive case, the stationary escape rate 
i st {x) is 



i s t(x) 



[a{p 3t {x))] 



1-H(x) 



r K x ) 



+ a(p st {x)) 



Pst(x). (40) 



One can see that the first term on the RHS of (|40[) is 
dominant for small ar and p.(x) < 1. For the jump 
density (HJ), in the limit a — > we obtain the stationary 
nonlinear Fokkcr-Planck equation 



d_ 

O.r 



dll 
2p—D{p at )p at {x) 



d_ 2 

dx 



mPW^)!. ( 4l ) 



where D (p st (x)) is the nonlinear diffusion coefficient de- 
fined as 



D(p st (x)) = 



a 2 [a( Pst (x))] l -^ x) 



•It, 



H(x) 



If we assume a zero flux condition at the boundaries of 
the interval [0, L], then 



dU 

J = 2/3—D{ Pst )p st (x) 



0_ 
dx 



[D(p st )p st (x)}=0, (42) 



and the stationary profile p s t{x) can be found from the 
nonlinear equation 



p st (x) = N- l D- 1 (p st (x)) exp {-2l3U(x)} 



(43) 



The time evolution of the density profile p(x, t) for the 
nonlinear fractional equation (|25p can be described as 

follows. At lower values of $ = J Q a(p(x, s))ds, the early 
evolution is the development of a single peak at the point 
of the minimum of p(x). This can be considered as an 
intermediate anomalous aggregation of particles. How- 
ever, incorporating the escape rate a (p) and the non- 
linear tempering factor e - * provide a regularization of 
anomalous aggregation. For sufficiently large $ the den- 
sity profile p (x, i) must converge to a stationary solution 
(|4"3l as t -> oo. 



IV. CONCLUSIONS. 

The aim in this paper was to derive the macroscopic 
nonlinear subdiffusive fractional equations for the evolu- 
tion of a mean density of random walkers by incorporat- 
ing a nonlinear escape rate and nonlinear jump distribu- 
tions. The main motivation was to take into account the 
interaction between particles on the mesoscopic level at 
which the random walker characteristics depend on the 
mean field density of particles. We showed that the non- 
linear escape rate leads to the effective regularization of 
standard subdiffusive fractional equations. Our modified 
fractional equations describe the transition from an inter- 
mediate subdiffusive regime to an asymptotically normal 
advection-diffusion transport regime. We showed that 



this transition is governed by a nonlinear tempering fac- 
tor that generalize the standard linear tempering. We 
discussed the aggregation phenomenon and showed the 
impact of a nonuniform distribution of anomalous expo- 
nent on aggregation patterns. 
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